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The temperature and field dependence of reversible magnetization have been measured on a 
YBa2Cu307_5 single crystal at six different doping concentrations. It is found that the data above 
2 T can be described by the scaling law based on the GL-LLL (lowest Landau level approach based on 
Ginzburg-Landau theory) critical fluctuation theory yielding the values of the slope of upper critical 
field — d_ffc2(r)/dr near Tc- This set of values is self-consistent with that obtained in doing the 
universal scaling for the six samples. Based on a simple Ginzburg-Landau approach, we determined 
the doping dependence of the coherence length ^ which behaves in a similar way as that determined 
from ^ = hvp/Esc with Esc the superconducting energy scale. Our results may suggest a growing 
coherence length towards more underdoping. 

PACS numbers: 74.40.-|-k, 74.25.Ha, 74.72.Bk, 74.25.0p 



In hole doped high temperature superconductors the 
transition temperature Tc and the maximum quasiparti- 
cle gap (or called as the pseudogap) behave in an oppo- 
site way: the former drops down but the latter rises up 
towards more underdoping 0. Although consensus has 
been reached on the doping dependence of some quan- 
tities, such as the transition temperature Tc, the super- 
fluid density ps and the condensation energy, etc., it re- 
mains still highly controversial about the doping depen- 
dence of the upper critical field or the coherence length 
in the underdoped region. In practice, however, to di- 
rectly determine i?c2(0) has turned out to be a difficult 
task due to its very large values. An alternatively way 
to derive — diJc2(r)/dT near Tc is to measure the re- 
versible magnetization or conductivity and then analyze 
the data based on the critical fluctuation theory. Using 
the Lawrence-Doniach model for layered structure of su- 
perconductors, UUah-Dorsey obtained expressions for the 
scaling functions of various thermodynamic and trans- 
port quantities around Tc Moreover, Tesanovic et al. 
pointed out that the scaling of magnetization due to crit- 
ical fluctuations near Hc2{T) can be represented in terms 
of the Ginzburg-Landau (GL) mean field theory on a de- 
generate manifold spanned by the lowest Landau level 
(LLL) i^j. By using a nonperturbative approach to the 
Ginzburg-Landau free energy functional, M(T) curves 
are evaluated explicitly for quasi-2D superconductors in 
a close form as: 
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where A and B are independent of T and H, but A is de- 
pendent on both the GL parameter k and \dHc2/dT\T^ , B 
depends on k. This scaling behavior is expected specially 
in a high magnetic field. Many experiments were tried 
to test these scaling laws and obtain the values of the 
mean- field transition temperature Tc{H) and the slope 
— difc2/dT ^, IMI^' However, due to the sample diver- 
sity, the scaling produced values oiTc{H) and — d//c2/dT 
that did not agree with each other. As a consequence, 
the universal scaling for superconducting diamagnetiza- 
tion fluctuations is still elusive. In this paper, we system- 
atically investigate the diamagnetization fluctuations in 
the vicinity of Tc of a YBa2 Cu^ O7-5 single crystal with 
six oxygen doping levels. The results indicate the plausi- 
bility of the existence of a universal scaling in the frame- 
work of 2D GL-LLL approximation theory. The doping 
dependence of i?c2(0) are also reliably obtained. 

The YBa2Cu307_5 single crystal used here was grown 
by top-seeded solution-growth using the Ba3-Cu5-0 sol- 
vent. Details for crystal growth were presented elsewhere 
7] . It has a shape of platelet with lateral dimensions of 
3.30 mm x 1.98 mm, thickness of 0.44 mm, and a mass 
around 15.56 mg. The different concentrations of oxygen 
were achieved by post-annealing the sample at different 
temperatures in flowing gas followed by a quenching in 
liquid nitrogen. The detailed annealing procedures are as 
follows: the as-grown YBa2Cu307_5 single crystal was 
first annealed at 400 °C for 180 hours with flowing oxy- 
gen and slowly cooled down to room temperature. The 
resulted crystal (SI) is close to optimally doped with a 
Tc = 92.0 K. Then the following doping status on this 
specific sample were achieved by annealing it in flowing 
oxygen in sequence: Tc — 85.3 K (520 °C for about 110 
h, S2), 79.5 K (540 °C for 120 h, S3), 68.5 K (580 °C 
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Table I. Parameters at six different annealed states 



sample 


Tc{K) 


P 




dHo2 
dT 


{T/K) 


SI 


92.0 


0.146 


92.7 ±0.6 


3.45 ±0.01 


S2 


85.3 


0.127 


85.7 ±0.3 


3.23 ±0.01 


S3 


79.5 


0.117 


80.0 ±0.5 


3.03 ±0.01 


S4 


68.5 


0.103 


70.1 ±0.5 


2.00 ±0.01 


S5 


58.6 


0.093 


62.6 ±0.5 


1.82 ±0.02 


S6 


30.5 


0.070 


36.7 ± 1.0 


1.43 ±0.05 



for 110 h, S4), 58.6 K (680 °C for 130 h, S5). The last 
sample (S6) was annealed at 520 °C for 90 h with flowing 
Na gas yielding = 30.5 K. 

The magnetization was measured by a Quantum 
Design superconducting quantum interference device 
(SQUID) magnetometer in both so-called zero-field- 
cooled (ZFC) and field-cooled (FC) modes with fields 
ranging from 10 G to 5 T parallel to c-axis. In the re- 
versible regime, the data measured using ZFC and FC 
modes coincide very well. In SQUID measurement, a 
scanning length of 3 cm was taken, the SQUID response 
curves in reversible regime were fully symmetric to avoid 
artificial signal. 

The T(.s' of the sample at six annealed states were de- 
termined by the deviating point of magnetization from 
the normal state background 8] (as shown in Fig.l (a), 
ZFC mode, 11=10 Oe). One can see that the super- 
conducting transitions are rather sharp near the tran- 
sition point. The gradually rounded foot of M{T) curves 
in more underdoped samples may be attributed to the 
easy flux motion since the system becomes more 2D like. 
The relationship between the TcS' and the oxygen doping 
level p [determined by using the phenomenological rela- 
tion Tc/r™'^'' = l-82.6(p- 0.16)2 taking T^'"^'^ = Qi.QK] 
has been summarized in Table I. 

Presented in Fig. 1(b) are the M(T) curves in the re- 
versible regime under the applied magnetic fields H =0.5, 
1, 2, 3, 4, 5 T. A crossing point at (T*, M*) appears 
in each set of M(T) curves. Such crossing behavior of 
the M{T) curves in high magnetic field has been well 
described by 2D or 3D GL-LLL scaling theory, and is a 
general consequence of fluctuations in the vortex state. It 
is interesting to note, however, that the value of M* /T* 
is not a constant in our six sets of M{T) data, be- 
ing contradicting with the theoretical prediction^. This 
deviation from the prediction of the 2D GL-LLL scaling 
theory were widely observed in Y123 [13, Hgl223 [ij 
and Bi2212 single crystals and may be attributed to 
the doping induced change of the anisotropy. 

Despite the deviation mentioned above, excellent 2D 
scaling curves are obtained for each set of data with 
different p and > 2 T, as shown in Fig. 2 where 
M{H,T)/(THf/^ is scaled as a function of the variable 
of [T - TXH)]/{THf/^. We also performed 3D scaling 
of M{T) for our six samples. The quality of the 2D scal- 
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FIG. 1: (color online)(a) Temperature dependence of the mag- 
netization measured in the warming-up process with a ZFC 
mode at = 10 Oe for the sample at six different doping 
levels, (b) Magnetization vs. temperature at magnetic fields 
of 0.5, 1, 2, 3, 4, 5 T. It is easy to see that a crossing point (or 
a narrow area) marked as (T*, M*) appears for each sample 
showing a typical behavior of critical fluctuation. 
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FIG. 2: (color online) Scaling curves based on the 2D GL-LLL 
critical fluctuation theory (see text) for the sample at six dif- 
ferent doping levels. From this scaling one can determine the 
slope \AHc2{T) / AT\t^ for each sample. Interestingly, a new 
type of crossing point exists for the scaled curves of different 
samples. 
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ing are, however, better than that of 3D scahng for the 
five sets of M{T) curves of S2, S3, S4, S5 and S6. For SI 
(close to optimally doped) the quality of 3D scaling is as 
good as that of 2D scaling with a narrower scaling region. 
This may suggest a 2D-3D crossover between p — 0.127 
and 0.146 for our deoxgenated YBCO crystals. 

To fulfill this 2D scaling, two variables Tco and 
dHc2/<iT are employed in Tc{H) as the fitting param- 
eters: Tc{H) = Tco - H{dHc2/dT)-^. The values of Tco 
and — di?c2/dT resulted from the fit are also listed in 
Table I. The critical fluctuation region ST = Tco — Tc in- 
creases with decreasing p, indicating a larger fluctuation 
regime for more underdoped YBCO. As shown in Table I, 
both the Tco and |diJc2/dr|7;- drops down towards more 
underdoping. Another interesting phenomenon shown by 
Fig. 2 for the scaled curves is that there is a new type of 
crossing point for different samples. This new type of 
crossing point may indicate a universal scaling among 
different samples. 

In the following we will check the feasibility of the uni- 
versal scaling for six samples. As all M{T) data are mea- 
sured on the same platelet of YBCO single crystal, this 
allows us to do universal scaling for six sets of M{T) 
curves. The vertical axis M / (TH)^-^ (y-axis) is intact 
because our experiments were done on the same sample 
with different doping concentrations and the GL param- 
eter K is weakly doping dependent. According to the 2D 
LLL-scaling theory, an analytical formula of the a;-axis 
for the scaling is written as 
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where L is a constant, related to s, the interlayer spac- 
ing. From this expression for x axis, it is evident that 
a full 2D scaling analysis of diamagnetization fluctua- 
tions has to include a material-dependent scaling factor 
A = ^ |^3t^It ■ -Based on the 2D LLL-scaling curve in 
Fig. 2, such full 2D scaling is performed by multiply- 
ing the "cc-axis" of each 2D LLL-scaled curve in Fig. 2 
by a factor A' = A{p)/A{p — 0.127) to make all curves 
collapse on a single branch. Here we use the values of 
obtained in 2D scaling analysis and leave k as 



dT 



a free fitting parameter. In doing the universal scaling 
we put the 2D-scaled data of 82 as a reference, i.e., the 
"a;-axis" is kept as it is. The resulted collapse of curves 
from the six samples is shown in Fig. 3. together with 
the doping dependence of k. It is found that the scaling 
quality is predominantly controlled by \dHc2/dT\T^ and 
K depends weakly on the doping level. Within the ex- 
perimental uncertainty the quality of the data collapse is 
reasonably good. For a comparison, a theoretical curve 
generated using Eq. (2) is also plotted in Fig. 3. This 
theoretical curve describes the data rather well except 
for a deviation for the sample p = 0.146 (Tc = 92.0) on 
which we believe the system becomes more 3D like. 
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FIG. 3: (color online)A universal scaling curve by collapsing 
the six scaling curves for different doping concentrations based 
on the closed-form of the scaling equation of Tesanovic et 
al.fl- The thick solid line is a theoretical curve after Eq.(2). 
The inset shows the relative k value derived from doing the 
universal scaling. 



We now discuss the results obtained by the universal 
scaling analysis. According to the Werthamer-Helfand- 
Hohenberg (WHH)[l3] theory for a 2D system, the value 
ofi/c2(0) isgivenbyi?c2(0) = 0.69TcX |di7c2/dT|T,. Dis- 
played in Fig. 4(a) is the dependence of i?c2(0) onp based 
on the obtained values of \dHc2/dT\T^ and Tc in Table I. 
Our data show a rough linear correlation between Hc2 (0) 
andp: i/c2(0) = 2620(p-0.058). We noted that such lin- 
ear -ffc2(0) — p was also obtained for underdoped Bi-2212 
polycrystal 14]. The spin-ordering quantum transition 
theory predicts that -ffc2(0) increases with doping in the 
underdoped regime 15]. Quantitatively, this theory cer- 
tainly deserves consideration in interpreting our data. 

It is interesting to note that in the underdoped regime 
the linearity of Hc2{0) — p leads to a linear depen- 
dence of superfluid density ps{0) on p by the relation: 
Hc2{0) = $o/27rC^ = $oK^/27rA2^. g^sed on the fact 
that K is weakly p dependent, which is indeed the case 
for our samples and also found previously in underdoped 
YBCO 1^, it is easy to have X^^iO) cx p,(0) cx p. This 
linear correlation was recently verified in Bi2212 [l^ . 
La2-2;Sr2;Cu04 and IIgBa2Cu04+5 [l^ . 

The coherence length ^afc(O) of each sample can also 
be extracted from the i?c2(0) value of Fig. 4(a). For 
example, we have Cab(O) = [$o/27ri?c2(0)]^/2 = 12.2 A 
for ffc2(0) = 220 T for the sample of Tc = 92.0 K. In Fig. 
4(b) £,ab is summarized and depicted as a function of p 
in the underdoped regime. As Hc2(0) is reduced to zero 
at the point p = 0.058 ± 0.002, very close to the critical 
point pc = 0.050 for superconductivity. This implies that 
i/c2(0) and Tc simultaneously drop to zero at the critical 
point of the phase diagram, indicative of the complete 
suppression of Ps(0) at pc ~ 0.05. Another consistent 
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FIG. 4: (color online) (a)The doping dependence of the up- 
per critical field -ffc2(0) (phase coherence) derived by doing 
the scaling. The solid line is a guide to the eye. (b)Doping 
dependence of the coherence length determined from two dif- 
ferent ways (see text). 



way to reckon the coherence length is using ^ — fivy /Esr 
with Esc — nksTc the superconducting energy scale jl8l|. 
vp is the Fermi velocity taking 2.5 x 10^cm/s[ljl| and is 
almost doping independent. In Fig. 4(b) we present also 
the coherence length calculated in this way with n « 20. 
It is clear that the coherence length obtained in these two 
ways coincide rather well. Our results here may indicate 
a growing coherence length in more underdoped region, 
being consistent with our earlier report [inj. 

In summary, we have systematically investigated the 
critical fluctuations on a YBa2Cu307_5 single crystal at 
six different doping concentrations. It is found that the 
data above 2T can be described by the universal scaling 
law based on the 2D GL-LLL critical fluctuation the- 
ory. Thus the values of the slope of upper critical field 
-diJc2(r)/dT near (and thus 7?c2(0)) can be reli- 
ably extracted. The coherence length derived from both 
Hc2{0) and ^ — Hvf/Esc show a similar growing behavior 
towards more underdoping. 
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